Local U(l) symmetry in Y(SO(5)) associated 
with Massless Thirring Model and its Bethe Ansatz 



Hong-Biao Zhang^'^ , Mo-Lin Ge^, Kang Xue^'^ 

1. Theoretical Physics Division, Nankai Institute of Mathematics, 
Nankai University, Tianjin 300071, P.R.China 
2. Educational Institute of Jilin Province, 
Changchun,Jilin,130022, P.R.China 
3. Physics Department, Northeast Normal University, 
Changchun,Jilin, 130024, P.R.China 



Abstract 

Keyword: Current algebra, U{1) gauge invariance, Yangian 
PACS number: 03.50.-d 03.65.Fd 



(I) Introduction 

Recently, it has been proposed by S.C.Zhang et al. that the antifcrromagnetic(AF) and 
superconducting(SC) phases of high-Tc cuprates arc unified by an approximated SO{5) 
symmetry principle[l]. Considerable support for this proposal came from numerical inves- 
tigations in models for high-Tc materials. In particular, it was shown that the low-energy 
excitations can be classified in terms of an S0{5) symmetry multiplet structurc[2,3]. Sub- 
sequently, extended Hubbard models and a two- leg ladder model related to 5*0(5) sym- 
metry have been introduced and analyzed in details[4,5,6]. On the other hand, Shelton 
and Senechal[7] have studied the problem of two coupled ID Tomonaga-Luttinger chains 
and concluded that approximate SO{5) symmetry can emerge in low-energy limit of this 
model. It is well-known that the Luttinger liquid is connected with the massless Thirring 
model. It is worth to deal with massless Thirring model with SO (5) symmetry. The model 
can be constructed by four-component fermions field operator ipi{x), we shall show that 
this model is exactly solvable by the Bethe ansatz method through a local U{1) transfor- 
mation, under which the fermion operator il^i^x) is transformed into q-deformed fermion 
operator ^i{x). This procedure leads to the diagonalization that looks in a simple man- 
ner by Wadati[8,9,10]. Furthermore, the study of Yangian algebra[ll,12,13,14] provides a 
significant tool in the formalism of integrable models. The generators of Yangian can be 
realized through currents for a given Lie algebra. It turns out that the current realization 
of 1^(5*0(5)) is not unique and allow a local U{1) gauge transformation. It is interesting 
to find the consequence of such a U (l)-freedom according to the q-deformation of fermion 
operator 

This paper is organized as follows: In the section (II), the massless Thirring model with 
SO{5) symmetry will be diagonalized and the Bethe ansatz wavefunction is constructed. 
In section (III), we shall give the current algebra realization of 1^(5*0(5)) in terms of 
q-dcformed fermion current that give rise to the local [/(l)-gauge transformation. 

(II) The massless Thirring model with SO (5) symmetry and its Bethe ansatz 
wavefuction 

Let us consider the massless Thirring model constructed by the four-component fermion 
field operator ip{x) = [ipi{x),'4>2{x),'4'3{^),'4'4{x)]'^ ■ The Hamiltonian takes the form: 



where Cij = Cji^Cu = and ni{x) = i/jf {x)'(pi{x) = 1,2,3,4) satisfy the anticommu- 
tation relations: 



For four-component fermion field operators tp{x) = [co-(x), d+(x)]^ and forms the cur- 
rent algebra obeying 50(5) [6]. In momentum space, this Hamiltonian can be written 
as: 




(1) 



[^+(x),V+(y)]+ = 
bPi{x),ijj{y)]+ = 
[tpi{x),ip^{y)]+ = dij5{x - y) 



(2) 
(3) 
(4) 
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+fyy j[ci2ct{k+i)c+{-k+^)c^{-k'+^)c^{k'+^) 

+Ci44(A; + |)4(-A: + l)d^{-k' + |)ct(A:' + |) 
+C23C+(fc + |)4(-A: + + + |) 



that obviously is made up of the pairs, so this model may be applied to superconducting. 
To diagonalize H, we introduce local U{1) transformation: 



+ C734d+(A: + h^ti-k + ^)di{-k' + ^)d^{k' + ^)]dkdk' dq (5) 



(^i{x) = exp[-i ^ 0ik(j)k{x)]^piix) (6) 
fe=i 

where 4'i{x) = j^^ijjf {y)'ijji{y)dy and Oik are constants. 

According to eq.(§)-eq.(§) and eq.{|6|) by direct calculation, we obtain(no summation 
over the repeated j): 

^i{x)^j{y) = -exp[ieij\^j{y)^i{x) (7) 
= -exp[i%]$+(y)$+(x) (8) 
= -ej;p[-i%]$+(?/)$^(2;) + 6ij5{x - y) (9) 
This is a special case of Zamolodchikov-Faddeev algebra[15,16]. 

en = {mod 27r) (10) 

6ijj- + Oji = (mod 2tt) (11) 

Therefore, eq.(p!o|) and eq.(|Tl|) are conditions given by the associativity of the special 
case of Zamolodchikov-Faddeev algebra. The meaning of eq.(|lO|) is clear that the particle 
itself must still be fermion for the same "i-spin" states, however, eq.(ll) show that the 
commutation relations between different "i-spin" states can be q-deformed and the q- 
deformation parameters should obey eq.(pTl) because of two-body interaction between 
different "i-spin" states. 

Under the local U{1) transformation eq.(^), the Hamiltonian eq.(|^) can be diagonalized 
and find the physical constraint conditions for real Cj and Cij. The Heisenberg equation 
idtiJi{x,t) = [ipi{x,t),H] reads: 

4 

dt'ipi{x,t) = vCidxipi{x,t) - i2g^Cijnj{x,t)'ilji{x,t) (12) 
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On account of the transformation eq.(|^ and the Heisenberg equation eq.(^) , we obtain: 

4 4 

dt^i{x,t) -vCidx^i{x,t) =i^[v{Ci - Cj)9ij - 2gCij]nj{x,t)ilji{x,t) exp[-i ^ Ojk(pk{x)] 

j=l k=l 

(13) 

By choosing 

% = ^^r^ ia^C,) (14) 

0^i = (15) 

then $j(x,t) satisfy the free-field equation . The Hamiltonian becomes diagonahzed (here 
we suppose: Cj ^ Cj; if Cj = Cj, the Hamiltonian can be diagonalized only when Cij = 0): 

^ r 

H' = ivJ2Ci / ^tix)d:,^i{x)dx (16) 
i=i •' 

The direct calculation shows that $j(x) and H' also satisfy the Heisenberg equation 
idt^i{x,t) = [^i{x,t), H'], so $i(x,i) are really dynamic variables regarding H'. 

Therefore, by using the local U{1) transformation eq.(|6|), the original Hamiltonian 
eq.(||) constructed by tpi{x) with anticommutation relations eq. (^)-eq. (^) has been trans- 
formed into the quadratic Hamiltonian eq.(|l^) in terms of the <5j(x) obeying q-deformed 
relations eq.(^)-eq.(|^. In the following, we shall show how the method in[8,9,10] works 
to find the Bethe ansatz wavefunction in a simple manner for SO {5) massless Thirring 
model. 

Let us denote by | ni,n2,n3,n4 > a eigenstate with rij <I>j-particles(i=l,2,3,4), it can 
be expressed by 

M ni 

ni,n2,n3,n4 >= 



...J Y[dxjip{xi,...,XM) n ^ti^jl) 



j=l jl=l 

n2 na n4 

X n ^^(a^Mi+ia) n HixAh+n) n ^ti^Ah+u) i o > (i?) 

i2=i i3=i i4=i 

where Mi = ni + n2 + ■■■ + Ui , M = and | > is the vacuum defined by 

iIj-[x)\{}>=Q (18) 

or equivalently 

^j{x) I >= (19) 
Substituting eq.(|l7|) and eq. ([l6|) into Schrodinger equation 

I ni,n2,n3,n4 > (20) 

it yields equation for 9j(xi, ...^xm) 

ivC^Ci )ip{xi,...,XM) = En^^n2,n3,niV>{xi,...,XM) (21) 

i=l j,=l ClXAh-i+h 
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whose solution is: 

M 

if{xi, ...,xm) = Aexp{i^kjXj) 

4 rii 

Eni,n2,n3,n4 = —v(^Ci ^ kM^.^+jf) (22) 
1=1 ji=l 

where kj and A are constants. Since the constant A is not essential, we shall omit it 
hereafter. The Bethe ansatz wavefuntion (^(xi, ...,xm) is defined by 



M 

ni,n2, 71-3,71.4 > 



„ „ M ni 

/ ... / Jl ...,xm) n 



i=i ii=i 

rj2 ns J14 

X [| Vj(a;Afi+i2) n '^ti^M2+j3) n V'|(a:M3+j4) I > (23) 

j2 = l i3 = l i4 = l 

Substituting eq.(^) into eg. (|T7|) , by detail calculation, we have 

I ni, 712,713, 714 >= / ... / Jl (ixj99(xi,...,a;M) n n n exp[i6lpg6'(j;Mp_i+jp-2;Af,_i+jJ] 

j=l l<P<q<4jp=li,=l 
ni 712 ns rii 

X n n ^ti^M^+n) n V'^(2;Af2+i3) n '^a^M^+j^) i o > 

ii=i i2=i j3=i i4=i 

A/ np Uq 

i=i i<p<g<4ip=iig=i 

ni n2 n3 n4 

X \[ il^ti^h) n V'J(a;Afi+i2) n V'3'(a;Af2+j3) n ^4'(a;A^3+i4) I > (24) 
ii=i i2=i i3=i i4=i 

where 0(x) = {if x < 0); 1 {if x > 0) and e{x) = 9{x) — 9{—x) hereafter . Thus, the 
Bethe ansatz wavefuntion (p{xi, ...,xm) takes the form: 

M n.p riq ^ 

(p{xi, xm) = exp[i J2 kjXj] Y[ n n [1 ~ ^*5^e(a;Afp_i+ip - XMq.i+jq)] (25) 

j=l l<P<q<4:jp = ljq = l 

which describes the many-body problem with (^-interactions. 

Suppose that M particles move in a region with the length L. For an arbitary 
Xj (Mp_i < J < Mp). Imposing the periodical boundary conditions (PBC), we have 



„ „ A/ rip nq Q 

J ... J Yldxj^p{xi,...,XM) n II I[['^-^i9^e{xMp.^+jp-XMq_^+jq)] 
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1 - itgepg/2 



kjL = —i 7ioZn -^^2- — h 2/o7r ih integer) (26) 



q=l 



I.e. 



kjL = — ^ riqOpq + 2/j7r {Ij integer) (27) 



q=l 
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that is exactly the Bethe ansatz equation. Obviously, the local U{1) transformation eq.(|6[) 
much helps the derivation of the Bethe ansatz condition for the massless Thirring model. 
(III). Current realization ofY{S0{5)) 

The SO{5) algebra does have the current realization, however the fermionic construc- 
tion is not unique. In parallel to the diagonlization of eq.(|l]) we shall show that the 
q-deformed operators shown in eq.(|6[) also provides a realization of S0{5) algebra, 

henceforth Yangian associated with S0{5). 

The original commutation relations of Y(g) were given by Drinfled[17,18] in the form: 

[Ix, Ifi] = Cx^ulu ; [Ix, J/j] = Cx^iyJu (28) 

[Jx, [Jfi, lu]] - [Ix, [Jii, Jiy]] = h^aXiiual3-f{Ia-, 1(3, I-y} (29) 
[[•^A) Jfj], [la, Jt]] + [[J(7, Jt], [h, J^l]] = h'^ {O'Xfiuaf^'yCaTU + acrTual3jCXf^u){Ia, 1(3, I-f} (30) 

where cxfj,u are structure constants of a simple Lie algebra g, h is a constant and 

OiXfiuafS-y — '^CXaaCfifSTCwypCcrTp i ■(^1;^2)^3} — ^ ^ XiXjXf^ (31) 

For Lie algebra SO{5), Y{SO{5)) is generated by antisymmetric generators {lab. Jab}- 
eq.(l^) reads 

[lab, led] = K^bJad + ^adhc " ^aJbd " ^bdhc) (32) 

[lab, Jed] = K^bcJad + ^adJbc — ^acJbd " hdJac) (33) 

lab = -ha, Jab = -Jba; (a, 6, c, d = 1, 2, 3, 4, 5) 

Not all of the relations in eq.(p9|)-eq.(pO|) are independent. After tedious calculation 
we can prove that there is only one independent ralation: 

fe, ^15] = ^/i'({/l3, ^42, 145} + {/l2, /45, /34} " [hi, ^42, h^} " {/l4, ^34, h^}) (34) 

where J23 and J15 are the Cartan subset. 

All the other relations other than eq.(p8[) can be generated on the basis of eq.(p^) by 
using Jacobi identities together with eq.( p2|) an d eq.(|3^). Therefore, for y(50(5)), eq.(p8D- 
eq.(^) can be expressed with eq.(^2|)-eq.(|34|) in such a simple manner. This conclusion 
can also be verified by RTT relation independently through tremendious computation. 

The generators of Y{SO{f))) can be realized by fermion current algebra as follows: 

lab = J Iab{x)dx ; Iab{x) = {x)T'''' -^{x) 

Jab = Tab + UJI^ ■ Tab = j (x)r''*5,^(x) 

Jab= / dxdye{x - y)Iac{x)Icb{y) (35) 



where F"* = — zF°F^, F*^ are 4x4 Dirac matrices, U = {h being arbitrary constant) 
and iIj{x) satisfies anticommutation relations eq. (0)-eq. (^) . It can be checked that the set 
{lab. Jab} Satisfies algebraic relations eq.(^)-eq.(p^ oiY{SO{b)). 
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As given by Ref[6] if V'(a^) = [ca{x), da{x)]'^ , then local generators Iab{x) of Lie algebra 
SO{b) are expressed in terms of spin S{x) = ^(c^(x)(Tc(x) + {x)(Td{x))^ charge Q{x) = 
\{c^{x)c{x) + d'^ [x)d{x) — 2) and 7f+(x) = —\c^ {x)aa2d^ {x) with 



Iab{x) 



\ 



( \ 

7rj^(x) +7ri(2;) 
7rJ(x) + 7r2(2;) -S^ix) 

+7r3(a;) 52(x) -Sx{x) 

Q(a;) i(7ri(x) - 7r^(x)) i(7r2(x) - 7rJ(x)) i(7r3(x) - 7r;^(x)) 0/ 

(36) 

where the value of matrix elements on the upper right triangle are determined by anti- 
symmetry, Iah{x) = -Iba{x). 

Under the local U{1) transformation eq.(^), the four-component fermion field oper- 
ators 'ip{x) = [ilJi{x),'tp2{x),il)z{x),'il)i{x)]^ is changed into q-deformed operator <I>(x) = 
Y^i{x)^(^2{x),^^{x),(^i{x)Y' ■ The generators of y(50(5)) is constructed by q-deformed 
fermionic current algebra as follows: 



lab 



Iab{x)dx ; Iab{x) 



Jab — Tab + U'J^ab 



T. 



ab 



^$+(x)r''$(x) 



ab^ J J dxdye{x - y)Iac{x)I cb{y) (37) 
Substituting eq.(37) into eq.(|32|)-eq.(|3^), we can obtain the constraint conditions: 

in - Ojn {mod 27r) (38) 



{i,j,m,n = 1,2,3,4) 

where eq. (|38|) sets the conditon making Y{SO{5)) constructed by q-deformed field operator 
$(x). This indicates that the current realization of 1^(50(5)) is not unique. Careful 
calculation shows that there are three free parameters in dij under conditions eq.(|l^), 
eq.(0) and eq.(38). So there exists an additional freedom in Y{SO{5)). Substituting 
eq.(U) into eq. (|3^) , we obtain 



Cirr 



a 



jn 



Ci Cn 



Ci Cn Cn Cn 



(mod 2tt) 



(39) 



{i,j,m,n = 1,2,3,4) 



In another words, H' is the Hamiltonian expressed by $i(x) , under the transformation 
eq(^, it becomes H where two-body interaction appears. So the physical meaning of U{1) 
transformation in Y{SO{5)) is connected with two body interaction in massless Thirring 
model with 50(5) symmetry. 

Prom the above analysis we see that there is a local U (1) gauge-invariance in the 
construction of the current algebra realization for Y{SO{5)). 
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Under local U{1) transformation eq.(^), eq.(36) is changed into 



Iah{x) 



( 

7f^'(x) +7fi(x) _0 

7r^(x) + 7r2(x) -Sj,[x) 

7f3"(x)_+ 7f3(x) 5*2 (x) 











(x)) i{'Kl{x) - T^^{x)) i{'K'i{x) - Txt,{x)) 0/ 

(40) 

where the value of matrix elements on the upper right triangle are determined by anti- 
symmetry, lahix) = -l^aix), and 



/ 5i(x) 
S2(x) 
V S^{x) 



I m^{^4>{x))zo^{^4>{x)) 
sin(^i±^0(x)) cos(^0(x)) 




-sin(^^±^(/.(x))cos(^(/.(x)) 
cos(^^±^0(x))cos(^0(x)) 
1 



Sx(x) 

S2(X) 

Sz(x) 



( 



+ 



■sin(^(/.(x))sin(^0(x)) 



cos(^^</.(x))sin(^^^ 




cos 



x))sin(^(/>(x)) 



c/>(x)) 



sin(4^(/.(x))sin(4^ 




0(x)) 




(41) 



W+(x) 
7r^(x) 



/ cos(^^0(x)) -sin(^^(/>(x)) 
2 ^ i sin(^i±^c/>(x)) cos(^<A(x)) 



vr+(x) 
7r+(x) 



(42) 



7f+(x)=exp(.^)[cos(^0(x))7r3+ + -A+(x)sin(^c/>(x))] (43) 

where 4){x) = X)i=i <Pi{x)i ^43 = /^i ^12 = a, ^'13 — ^42 = i^, SC order paremeter A'^{x) = 
—ic'^{x)a2d'^{x) and AF order parameter iV(x) = ^(c"'"(x)fTc(x) — d'^ {x)ad{x)). 

From the above analysis we see that under the U{1) transformation eq(^), eq.(^) 
still obey Yangian algebra as eq.(|35|) does if 6ij satisfy the conditions eq.(|To|),eq.([Tl|) and 
eq.(^). This indicates that there is a local ^7(1) gauge- invariance in the construction of 
the current realization for Y(S0{5)). It turns out that after the transformation, there are 
local phase factors in the current realization of Y{SO{5)) (shown by eq. (pl|) ,eq. (^2|) and 
eq.(l^), but eq.(|37|) still satisfy Y{SO{5)) Yangian algebraic relations, i.e., there is a local 
C/(l) gauge-invariance in such a current realization of Y{S0{5)). 

We also find that the transformation eq.(|6|) can be used to diagonalize the massless 
Thirring model with 5*0(5) symmetry that will help to understand the physical meaning 
of the introduced local U (1) symmetry. 

In another words, H' is the Hamiltonian expressed by $i(x) , under the transformation 
eq.(^, it becomes H where two body interaction appears. So the physical meaning of 
U{1) transformation in the current realization of 1^(50(5)) connected with two body 
interaction in this physical models. Applying the transformation we find the local U{1) 
gauge-invariance in Y{S0{5)) explicitly . 

Noting that there are some non-trivial phase factors in the generators of Yangian, 
but they still satisfy the commutation relations of Y(SO{5)), i.e., there is a local U{1) 
gauge-invariance in Y{S0{5)). 

(IV). Conclusion and Acknowledge 
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Using a local U{1) transformation connecting the four-component fermionic field op- 
erator Vi(^) with q-deformed one , it is helpful to diagonalize massless Thirring 
model with SO (5) symmetry. The Bethe ansatz wavefuntion is obtained in a simple man- 
ner. It turns out that the current realization of Y{SO{5)) is not unique and exist a local 
U{1) gauge transformation. This shows the existence of a local U{1) symmetry in the 
current realization of Y{SO{5)). Correspondingly, the transformation leads to the local 
C/(l)-gauge invariance for Y{SO{5)). The explicit forms of phase factors for 5*0(5) has 
been shown. 

The authors would like to thank Dr.Jing-Ling Chen helpful discussion . This work is 
in part supported by NSF of China. 
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